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^— s ' In the context of DifTerential Renormalization, using Constrained Differential Renormalization rules at 

one loop, we show how to obtain concrete results in two loop calculations without making use of Ward 
identities. In order to do that, we obtain a list of integrals with overlapping divergences compatible with 
Oh. CDR that can be applied to various two loop background field calculations. As an example, wc obtain 

(— ! , the two loop coefficient of the beta function of QED, SuperQED and Yang-Mills theory. 
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1 Introduction 

Since its introduction in 1991 [1] differential rcnormalization (DR) has been applied in a large variety of 
situations like (j>* at three loops [1], Yang Mills at one loop [1], QED at two loops [2], the Wess-Zumino 
model at three loops [3], SuperQED at two loops [5] and Super YM at two loops [4] among others. In 
all of them the method has exhibited one of its strengths: the simple and elegant way of obtaining the 
coefficients of the renormalization group equation (beta functions and anomalous dimensions). 

The classical objection to the method points out that the renormalization of each divergence introduces 
an "a priori" independent scale. In order to preserve gauge invariancc one has to apply at each order 
Ward identities in order to obtain the correct relations among scales. A constructive set of rules that 
avoids this somewhat tedious procedure has only been deviced at one loop. The enhanced set of rules that 
goes by the name of Constrained Differential Renormalization (CDR) [6, 7, 8, 9] automatically produces 
expressions with a single mass scale that fulfill Ward identities. 

At two loops a fully consistent extension of CDR is not known. However we believe that the results 
of this paper are a small step ahead in the right direction. Mainly, we will be dealing with the divergent 
parts of the amplitudes. For them we will show how one-loop CDR still has the strength to produce gauge 
invariant answers. Interestingly enough, from the divergent parts useful results can still be obtained such 
as the beta function and the anomalous dimensions. 

The structure of the paper is as follows. First we briefly review the basic rules of DR and CDR. In 
the next section we write a list of renormalized expressions for several two-loop integrals that contain 
overlapping divergences and comment about the general strategy to obtain them. As an example of the 
use of these integrals, we obtain the two loop coefficient of the beta function of QED, SuperQED and 
Yang-Mills. The first two calculations basically reproduce the results obtained in [2] and [3] in a more 
efficient way. 

2 DR and CDR 

2.1 Differential Renormalization (DR) 

DR is a renormalization method that consists in replacing coordinate-space expressions that are too 
singular by derivatives of less singular ones. This method does not need cutoff nor explicit counterterms, 
although they are implicitly used when performing formal integration by parts. The basic idea is that 
divergent one loop expressions are well defined for non-coincident points, but at short distances the 
amplitude is too singular and does not have a Fourier transform. So, to renormalizc one can replace 
the divergent expression for the derivative of a less singular one that has the same values as the original 
outside the origin, and has a well defined Fourier transform (if formal integration by parts is used with 
the derivatives). 

This method is especially well suited for supersymmetric theories because we stay in four dimensions 
all the time, which is not the case with other methods as dimensional regularization or reduction. 

As an example consider the one loop contribution of (f>* theory. The bare expression is 
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Both expressions coincide for x ^ 0, but the new one has a weU defined Fourier transform if we neglect 
the divergent surface terms that appears upon integrating by parts the d'Alembertian. Here is where the 
counterterms hide, and by applying formal integration by parts [1] we are implicitly taking them into 
account. I.e. 
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A constant with mass dimension M has been introduced for dimensional reasons. It parametrizes the 
local ambiguity 

Inx^AI'^ Inx^A'P M' ,, , 

A crutial observation is that this shift M -^ M' can be absorbed in a rescaling of the coupling constant A 
[1]. This is a hint that renormalized amplitudes satisfy renormalization group equations, with M playing 
the role of the renormalization group scale. 

Although we will not treat massive theories here, they have been considered in [1] and [10]. The 
appearance of a bare mass does not interfere with the method since the DR is related to short-distance 
singularities and masses only change the long-distance behaviour of the correlators. 

As was previously mentioned, each divergency is cured independently by a suitable replacement like 
the one above. An important issue arises with symmetric theories, where gauge invariance should be 
satisfied at each order in perturbation theory. The standard lore is that Ward identities (except for 
anomalies) can always be satisfied with the renormalized expression by adjusting the different mass 
scales obtained for different diagrams and expressions (ie. fixing a renormalization scheme); but at the 
same time, one always has to explicitly apply these identities to obtain the correct scale relation. 

2.2 Constrained Differential Renormalization (CDR) 

The basic idea is to obtain a renormalization scheme that automatically preserves Ward identities for 
renormalized amplitudes. For this, the following rules will be applied 

1. Differential reduction 

• Functions with singular behaviour worse than logarithmic are reduced to derivatives of (at 
most) logarithmically divergent functions without introducing extra dimensionful constants. 

• Logarithmically divergent expressions are written as derivatives of regular functions, introduc- 
ing one single constant M, which has dimensions of mass and plays the role of the renormal- 
ization group scale. 

2. Formal integration by parts We do not take care of the surface terms that appear when applying 
integration by parts. Related to this, differentiation and renormalization must be two commutative 
operations: let F an arbitrary function, then [dF]ii = d[F\R. 



3. Renormalization rule of the delta function: 

[F{x, xi,..., Xn)S{x - y)]R = [F{x, xi,..., Xn)]RS{x - y) (2.4) 

4. Validity of the propagator equation 

[F{x,xi,...,Xn)iO - m^)A„(x)]fl = -[F{x,xi,. . . ,x„)(5(a;)]_R (2.5) 

where A™ is the massive propagator of a particle of mass m and F an arbitrary function. 

Central to the fulfillment of the Ward identities (and the action principle, from which they can be 
derived) is that the application of the kinetic differential operator to some propagator line inside a 
Feynman graph is equivalent to the contraction of the line to a point. This statement is guaranteed 
to hold through the previous set of rules. The upshot is a basic set of renormalized expressions (basic 
functions) with different numbers of propagators and various differential operators acting only on one of 
the propagators, and involving a single scale M. Therefore the CDR program amounts to the following 
two step operation: 

• Express the Feynman diagram in terms of these basic functions performing all the index contractions 
(this is an important point, because CDR does not commute with index contraction) and by means 
of the Leibniz rule moving all the derivatives to act on one of the propagators. 

• Replace the basic functions with their renormalized version. 

Applying the previously defined rules, one can obtain a set of renormalized basic functions. Here we 
present a list of the most relevant ones (A stands for the masslcss propagator A = j^^)- 
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TR[d^d,] = TR[d^d,^-S^,n] + -S^,TR[D]~^^^^S^J{x)Siy) (2.7) 



Out from the list of CDR results for three propagators [6, 8], in our calculations we will make use of 
the decomposition into trace and traceless part that is imposed by CDR as 

A ^ J I ^ A"> KL J -^287r2 

where we have defined T[0] = AAOA. 

When using other gauges different from Feynman gauge, some bare expressions are written in terms 
of a quantity we define as A{x) = ^/^^'is Ins^g^^ where s is an irrelevant constant with mass dimension. 
For this structure, CDR prescribes [6] 

{AuA)r = A| 
{Ad,d.A)n = i(^A?,,-^a,9.ij) (2.8) 

CDR has been checked in abelian and non-abelian gauge theories [6, 9], and in supersymmetric 
calculations [11, 12]. 



3 Two loop use of one loop CDR results 

After CDR has been applied, all the ambiguities (local terms) at one loop are fixed in such a way that 
Ward identities are automatically fulfilled; we will now see how well we can do in two loop calculations. 
We anticipate that at two loops we will be only interested in the divergent part. 

Before we proceed with all of this, let us set up some notation that we use in this paper. In our 
calculations we have found that most of the results can be put in terms of an integral expression that we 
designate as I^ 

iHx - y) = J d\A,^Al^ (3.1) 

which is renormalized applying CDR rules as 



1 Inx^RP 
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4(-) = tt:^--^'^ (3-2) 



3.1 Nested divergences 



This case is particularly simple because CDR can be applied in a systematic way. Starting from the 
"inner" divergency, its regularization according to CDR gives an expression with logarithms of a single 
scale (Inx^M^) and fixed local terms. The one-loop Ward identities are fulfilled. In the next step, when 
tackling the outer part of the diagram, a simple logarithm like the one shown above is promoted to an 
expression of the form In^ x^M^ + m In x^AP, with m a calculable coefficient; at the same time, the finite 
terms multiply outer divergences which will produce additional logarithms of new scales (the genuine two- 
loop scales). CDR does note yet prescribe what the different scales of this later step should be. Hence, 
we may take all of them the same, and equal to M, at the price of leaving undetermined finite terms. In 
summary, as long as we are interested in two-loop quantities that just depend on the logarithms, we see 
that CDR at one loop prescribes a unique answer. If we were interested in three-loop calculations, the 
two loop finite parts would be essential. 

This simple scheme has some subtleties when studying diagrams with indices, because index contrac- 
tion does not commute with CDR. Therefore, the correct order is to insert into the outer diagram the 
non-renormalized expression for the "inner" one loop diagram, perform all the index contractions, and 
then renormalize. 

With this procedure we have renormalized all the different structures made up with I^ that we have 
found in our calculations. The divergent parts of those structures are (with . . . standing for the finite 
local contribution that we are not taking into account) 
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3.2 Overlapping divergences 

In this section we will focus on the subtle point of the application of CDR results to a two loop calculation: 
diagrams with overlapping divergences. In these cases it is sometimes difficult to recognize the one loop 
subdivergences that need to be treated with CDR to start with. The idea will be to obtain through 
different methods a list of renormalized two loop integrals with overlapping divergences where in each 
calculation one loop CDR rules have been maintained in every step. 

The list is restricted to integrals with at most four derivatives acting on the propagators and two free 
indices. This is a basis with which we can express all the overlapping two loop contribution to two point 
functions in most theories with derivative couplings in a background field approach (as Yang-Mills, QED, 
SuperQED or SYM). 

We use the conventions oi z — x — y and 9^ = (9j^ . We also define H as 

H[0i,02 ; 03,0i] ^ j d\d\ (Of A,,)(O^A„,)(0|Aj,„)(0|A,„)A„, (3.4) 

being Oi a differential operator. Finally, let us remark that as with nested divergences, we will not 
obtain the finite parts of these expressions (we are only interested in the divergent parts, which are 
unambiguously fixed by one loop CDR). In the final results . . . will stand for these terms. 
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iln2z2M2 + |lnz2A/2 



(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 
(3.11) 

(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 

(3.18) 
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(3.20) 



These integrals are obtained basically by applying two properties: 

• Integral relations presented in appendix D. These exact relations allow us to put some of the 
integrals in terms of others that have an explicit d'Alembertian acting on one of the propagators. 
Once we have done that, we can apply DA = —S to put these integrals in terms of the previously 
defined I^. Then, we can straightforwardly apply the procedure for renormalization with nested 
divergences that preserves one loop CDR.^ 

• The decomposition into trace part, traceless part and fixed local term imposed by CDR to T[d^d,y] 
as (2.7). 

Let us show in an explicit example how we apply these two procedures. Considering integral (3.9), 
this can be evaluated with both methods. First, we will make use of (D.2) and put this integral as sum 
of different integrals that have the divergences nested 
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/ d^ud^v A^u^xv{dld^Ayu)AyyAuy = 
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+ /<i*"<i'.A...A,.pi.,)Wi,.)i..- 
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(3.21) 



Now, we have to apply DA = —b and rewrite this integrals in terms of I^ (note that a = -j^^t^\ 
The third integral can be easily shown to be finite, and its value is obtained in appendix A. 
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\ d^ud'^v AxuA^^{d\dlAyu)Ay^A^ 



-ia.(A/i) + i(Aa./i) + ^9,(nA) 
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Applying the renormalization procedure for nested divergences, the final renormalized expression is 
found to be 



dll d^udh, A,^AUdldyAyu)Ay,A, 



d,,n^ :r-^ + ... (3.23) 
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Where . . . stand again for the finite terms that we are not taking into account and z = x — y. 

The other method of obtaining this integral is making use of the CDR relation (2.7) and perform a 



^Also this is the reason why we have not Usted here the cases where the differential operator is a d'Alembertian. For 
example, it is obvious that -ff[D, 1 ; 1, 1] = — A/^. 



trace-traceless decomposition of (9^9^Ayti)Ay„A„„ adding up the CDR term, i.e. 
Off d^ud^v A,„A,„(a^9^Ay„)Aj,„A„„ = 
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(3.24) 
(3.25) 
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Where d\Ix,y b. is the traceless part, that is finite. As we can see, both results agree. 

Although in this example we can perform the calculation with both methods with the same effort, 
with other integrals the situation is different, being necessary to study each expression in order to choose 
the best one. The explicit evaluation of all the integrals is presented in appendix A. 



4 Abelian examples 

The two loop renormalization of QED was carried out in [2], whereas SuperQED was studied in [5]. 
We will reobtain the two loop beta function coefficient of those theories using the procedure previously 
sketched, so that no Ward identity will be used. 

We will make the calculations in both cases in a background field approach [15], which implies that 
we will only treat the background field selfenergy (it is the only function needed to be considered when 
obtaining the beta function [15]). 



4.1 QED 

In this calculation, we apply the definitions and conventions of [2] . The only difference is the background 
splitting in the gauge field of the form A^ -^ A^j + i?^ (similar to the Yang-Mills case [15]). 

4.1.1 One loop 

For completeness, we briefly review the CDR renormalization of the photon selfenergy 11^,^ [6]. Although 
we need this amplitude to obtain the one loop coefficient of the beta function, it will be of no use when 
obtaining the two loop one. The bare vacuum polarization is 



Tt(1 loop) 



^{lefTr [7^7^93^A7,7-'a|'A] . (4.1) 

Simplifying the Dirac matrices and applying afterwards CDR we arrive to the renormalized expression 
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Figure 1: Two loop QED diagrams 



4.1.2 Two loop 

Proceeding with the two loop case, there arc two relevant graphs with external background fields. Diagram 
(a) has the divergences nested, whereas diagram (b) has overlapping divergences (figure 1). 

• Diagram (a) 

The expression for this diagram is 



^t'-Kx) = -(* 



d ud V Tr 



7^7'(-ajA,„)s(i)(u - «)7^(-a,"A„,)7.7'^(-5^A,,)' 



(4.3) 



where E^^^ is the one loop fermion selfenergy. As the term that takes care of the running of the 
coupling constant in the RG equation is not relevant when obtaining the two loop beta function [2], we 
could restrict ourselves to Feynman gauge. In that gauge, the bare fermion selfenergy is 



Y.^^\x) = -2eVA9AA(a;) 



(4.4) 



As we stated in the previous section, CDR imposes a strict order to the operations of index contraction 
and renormalization: First all the indices should be contracted, and only after that we can renormalize. 
Inserting the bare fermion selfenergy we are keeping this order here. 

Applying the Clifford algebra of the 7-matrices and the integral indentity 



/ d^ud^v AxuAyy{AuvdxAuv) = -^df / d'^ud'^v A^uAy-uA. 
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uv 



is straightforward to write 11^,^°' in terms of I^ as 



n(2/)(a;) = 8e4(9^A)(a,/i)-4eV(5AA)(9A/') 



(4.5) 



(4.6) 



We only have to replace the different structures made up with I^ by their renormalizcd values as were 
listed in section 3.1. The final renormalized expression is then 



n?.1(^) 
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Sf^^an 



In^x^il/^ + flnx^M^ 



(4.7) 



• Diagram (b) 



This diagram, opposite of the previous one, has overlapping divergences. We will express this in terms 
of the listed overlapping integrals. The bare expression is 

^%'\^ -y) = -{lef j d^ud^Tr [7^(-7"aSA,„)7''(~7^a;^A„,)7,(-7^9];A,„)7p(-7'^9.A„,)A„,„] 

(4.8) 

Applying algebra of Dirac matrices we can expand this expression, the complete result being listed in 
appendix B. As can be seen there, the contributions with the delta (DA = —5) can be easily expressed 
in terms of /^, and its renormalization is as in the previous case. The other contributions can be found 
in the list of renormalized expressions with overlapping divergences (or can be easily expressed in terms 
of integrals of the list). So that, we only have to apply directly the results found in section 3. No Ward 
identity is needed to fix the correct scale. We found the divergent part of the renormalized expression to 
be 

nl?](.) = 



12(47r2)3 
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-\n'x'M'-flnx^M' 



■S^^oa 



In'^x^A'P + ^lnx'^AP 



(4.9) 



• Final expression 

With the two previous results, the total two loop renormalized background selfenergy is 
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Inx'^AP 



(4.10) 



where again . . . stands for the local terms that we are not taking into account. We can use this result 
(along with the one loop renormalized value (4.2)) to obtain straightforwardly as in [2] the two loop 
expansion of the beta function in terms of a = |:^ as 



/3(«) 



+ -.+0(a^) 



(4.11) 



2q2 

~3^^ 27r2 

To stress the key points of our calculation, let us compare this procedure with usual differential 
renormalization [2] . Being Ms and My the one loop renormalization scales of the fermion selfenergy and 
the three point vertex V^ respectively, the Ward identity 



d r 1 

— ■V^(.T- z,y-z) = -ie S'-'^^z - x) - S'-^^z - y) T.{x - y) 



imposes that these scales are related as [2] 



, Ml 

In — If 

Ml 



(4.12) 



(4.13) 



When dealing with the two loop contributions, in each case we have to make use of the corresponding 
one loop scale {My: or My), being found the final values for HI,,^' and Tllii, to be 

Inx^M^ 



Ii?:^\{x) = - 



-^/ji/ R\ 



n 
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fiv R 
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967r2 Vtt 
1 /a 
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In^ x^M'l 



I In x2 A/2 



- V°°- 



'inx^AP + ^lnx^M^^ 



-SnuOa 



Inx^M^ 



(4.14) 
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Then, to obtain the enthe two loop vacuum polarization we have to use the mass relation (4.13) to 
put one of the scales in terms of the other. As can be seen, all of this is avoided in our procedure. No 
Ward identity was needed to be imposed, because CDR at one loop has fixed all the ambiguity relevant 
to logarithms of the scale at the two loop level. 

4.2 SuperQED 

We will now study another example of the use of one loop CDR results in a two loop calculation. In 
this case, we will discuss the supersymmctric extension of the previous case, SuperQED. This example 
was already studied in the context of differential renormalization in [5], although as in the previous case. 
Ward identities were imposed to obtain the correct scale relations. In this calculation, we will use the 
supersymmctric conventions for the TV = 1 superspace of [16]. With them, the propagator is formed by a 
product of the delta function of Grassmanian variables and the usual part Pij = 5{6i — 6j)l^{xi — Xj). The 
action, expressed in terms of a chiral abelian superfield Wa = iD^DaV (where Da is the supercovariant 
derivative and V a real superfield ) and chiral matter superfields <&+ and <&_ is of the form 

S = f d^xd^ew^+ f d^xd'^e^+e^^"^++ f d^xd*e^-e^^<i>- (4.15) 

In this case the background splitting is defined a,s V ^ V + B [5] . 

A complete study of perturbative calculations in superspace can be found in [16]. The basics of the 
method is to express the diagram in terms of superpropagators and supercovariant derivatives acting on 
them. Integrating by parts, we have to make the derivatives act on some propagators and afterwards 
apply the property 612D1D1612 = <^i2- Here we will write directly the expression found after this step, 
being the previous standard superspace calculation found in [5] and [16] for the SQED case. 

4.2.1 One loop 

For completeness, as in the QED case, we consider the one loop renormalization, although we will not 
use it in the two loop calculation. The bare amplitude is 

r("°°P) = T / d''xd^yd''OB{x, 0) {D''D^DaB{y, 9)) A^^ + 

+ ^ f d''xd^yd''9Bix,e)Biy,e)A,ynA,y (4.16) 

and renormalizing according to CDR rules we found 



r^' '°°^' = "16(4^ / '^"^^'2/'^"^ ^(^' ^) {D"D'DaB{y, 6)) D- 



{x-yy 



(4.17) 



4.2.2 Two loops 



We first list the bare expressions of diagrams presented in figure 2. We will write these expressions in 
terms of I^ and integral (3.13). Then, we renormalize them according to the procedure defined in section 
3. 

• Diagram (a) 
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Figure 2: Two loop SQED diagrams 



p(2 loop) 



• Diagram (b) 



-,(2 loop) 



• Diagram (c) 



f d'^xd'^yd^e B{x,e)B{y,e) [□(A/1) - 2A3 - a""(Aac,d/^)] 



+ y / d^xd^yd'^e B{x, 9) {D°'D^Da.B{y, 6)) A/^ 



g* I d'^xd'^yd'^OBix, 6)B{y, 9) [-□(A/^) + 2A3 + a""(A9„^/i)] 
/ / d^xd'^yd'^e B{x, 9) {D^D^D^Biy, 9)) AI^ 



(4.18) 



(4.19) 



^i' '°°P) = Y d^xd^yd^9B{x, 9) {D''D^D^B{y, 9)) [AI^] + 

4 r 

+ ^ d'^xd'^yd'^9 B{x, 9)B{y, 9) [□(A/^) - A''' - 8°"^ [Ad^o.!^)] + 

4 f 

+ ^ d*xd*yd'^9 B{x, 9) {D'^D^D"'B{y, 9)) C"^" H[9^^, 1 ; 1, a„^] 



Where C"'^ is the raising and lowering matrix of the SU(2) supersymmetric indices. 
• Diagram (d) 



n(2 loop) 



4 p 

-^ / d^xd^yd'^9 B{x,9)B{y,9)A^ 



(4.20) 



(4.21) 
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• Final renormalized expression 

We can now substitute tlie previous expressions by its renormalized values and add tlicni in the final 
result. It can be easily shown that only the the last part of (4.20) survives, and the rest cancels exactly. 
But we are allowed to do this only because we are using CDR at one loop. In [5] this simplified calculation 
was forbidden. One has to renormalize in a separate form each diagram; obtain the renormalized results 
with different scales and relate them via Ward identities. Only at this point we can add up the expressions. 
With our method, we always renormalize the same structures in the same way and with the same scale, 
so that we can add up the expressions even before performing the renormalization. 

The final result is then (multiplying by 2 because we have performed the calculations with only one 
of the chiral matter fields $+, $_). 

^(2 loop) _ „ /^{2) p(2) p(3) p(4) \ 

^i?, — "^ y'^ a R + '^ b R + '^ c R + '^ d Rj 

= 9^ f d'^xd^yd'^e B{x, 9) [D^D^D°'B{y, 6)] C^^ i?^[9;3^, 1 ; 1, 9aa] 



16(4.2)3 / d'xdSd'9B{x, 0) {D'^D^D^Biv, 9)) U ^^ _ ^^^ 



(4.22) 



Where we have used the result (3.13) of the list of integrals with overlapping divergences. Taking 
into account that with our conventions (those of reference [16]) g — \igsQED; with qsqed the usual 
coupling constant of SupcrQED, we can easily evaluate the two loop expansion of the beta function as 
in [5] to be 

f3{gSQED) = ^9SQED + ^^9SQED + OiglgEo) (4-23) 

Which agrees with previous results found in the literature [13, 14]. 

5 Non abelian example 

One of the relevant calculations that was not yet obtained with differential renormalization is the two 
loop renormalization of Yang-Mills theory. Here we will detail how we can obtain the two loop coefficient 
of the beta function of this theory with little effort, making use of the procedure defined in section 3. 

We will follow [15] where a two loop calculation of the beta function of Yang-Mills theory in the 
background field approach with dimensional renormalization was performed. The Yang- Mills lagrangian, 
written in terms of gauge covariant derivatives {D^" = df,S'"'+gf'''"'B''^ and D^^^ = d^d'"'+gP'"'{B^+A''^)) 
and ghost fields rj"- is 

^ = \F;^F^, + ^{D,A,r{D,A,r + {D,f)nB,vr. (5.1) 

We have to point out that the gauge fixing parameter a will be redefined in our calculation as — = 1+^^, 
so that Feynman gauge (a = 1) will correspond to ^ = 0. 



5.1 One loop renormalization 

We now briefly review the results found in [9] and in [1] for the one loop differential renormalization of the 
c^uantum and background fields. Although the background selfenergy is all that we need when obtaining 
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the one loop beta function, we also have to consider the renormalization of the quantum selfenergy as 
this will be a one loop insertion in a two loop diagram (diagram (b)). Also, it will be needed in order 
to obtain the first coefficient of the expansion of the function that takes into account the running of the 
gauge parameter in the renormalization group equation. We will detail this later in section 5.3. 

5.1.1 Correction to the 5° propagator 
(Feynman gauge ^ = 0) 

This is the sum of two different diagrams: one with a loop of quantum gauge fields, and another with 
ghost fields (see figure 3). 

The procedure is as defined in the CDR section. First we write the expressions in terms of the basic 
functions, and after that we replace them by their renormalized values. 

• Gauge loop 

Here and in the rest of the following diagrams of Yang-Mills theory, Df^^ denotes a space-time deriva- 
tive acting on one external field. Applying the Leibniz rule, this becomes a minus derivative acting on 
the propagators. 

^2 racd J'bdc 



< B;ix)Bt{y) > 



2 ^"^^ 



d,A^ - 8S^,nA^ + 8df,iAd,A) - 16A9^9,A] 



-26,pDy + 6p,{dy- dl) + 2b,,Dl 



(5.2) 



• Ghost loop 



< Bl{x)Bl{v) > = -g'r'\f^''A^y{dl -dl){dl- dl)A,y 
= -g^CA5''''[2dp{Ad,A)-AAd^d,A] 



(5.3) 



• Total contribution 

The total non-renormalized contribution is 

aab 



< Bl{x)Bl{y) > = g^CAb'^'' [^d^d.A^ - AS^^UA^ + 2a^(Aa,A) - 4A9^a,A] 



(5.4) 



Replacing the values prescribed by CDR, the renormalized one loop contribution to the S" propagator 
is obtained as 



<B'^(x)Bl{Q)>R = g^CA5''\d^,d,-5f,,n) 



g^CAS^'id^d, - S^^D) 



11 a2 1 c. n 

—Al ttS(x) 

3 ^ 727r2 "■ ' 



11 Inx^M^ 

-D- 



487r2(47r2) x^ 



727r: 



Ax) 



(5.5) 



The result found here is transverse, fulfilling the corresponding Ward identity as was guaranteed by 
the use of CDR. 
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vAA> 



sSSL 



Figure 3: One loop YM diagrams 

5.1.2 Correction to the A" propagator 
(Feynman gauge ^ = 0) 

As in the previous case, we have contributions with gauge and ghost loops. Proceeding in the same way 
as in the background selfenergy, we first obtain the full expanded bare expressions and then renormalize 
them according to CDR rules. 

• Gauge loop 



2 facd £bdc 



< Alix)AUy) > 



gljacaj 



-A, 



5^,{Dl- dl) + 6.,{d; - d;) + 5p^{dl ^dl) 



s.AD^p - d'ji) + Sp.idy ~Di) + 5,„{dy~ di) 



g^CAS"'' 



[2(8^8, - ^^.□)A2 + 109^(A9.A) - lOAa^a.A- 



iSf,,d\A8xA)-2S^,AinA)] 



Ghost loop 



(5.6) 



<Alix)Aliy)> 



-9 



2 j^adc j^bcd 



f^'^A^y 81 81 A 



xy 



= -g^CAS'''[d^{Ad,A)-Ad^d,A] 



(5.7) 



• Total contribution 

Adding the two previous results we found the non-renormalized contribution to be 

<A;{x)Al{y)> = g^CAS''''[df,d,A^-6^,nA^ + 4df,{A8,A)-2S^,d\A8xA)- 
- 4Aa^a,A - 5^.A(nA)] 

and with CDR identities it is straightforward to obtain the renormalized contribution as 



<A';,ix)AtiO)>B = q'CaS'- 

^gHJAS^ 
1447r2 



1(8^8, - 5^u^)Al - ^{8^8, - S,^a)6{x) 



{d^,8y - 5fj,yU) 



727r2^ ^ 
15 hix'^M'^ 
47r2 ^2 



25{x) 



(5.8) 



(5.9) 
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5.1.3 Effective action in a generic gauge 



All of this calculation is performed in Feynnian gauge, so that we have to take care of the running of the 
gauge parameter. In order to do that, we consider a functional approach. 

The idea is to expand the effective action at one loop at second order in the background fields and 
retain only the linear dependence in the gauge parameter, due to the fact that in the renormalization 
group equation we will derive wrt. this parameter and after that impose Feynman gauge (^=0). 

To obtain the exact one loop effective action it is well known that we have to consider only the part 
of the lagrangian quadratic in the ^" fields. This part is 



where □"'' = [Df'D^Y^. 

Then, the generating functional for connected Green functions can be put as 

1 



W 



-trXn [J^.n-" - 2gr'Bl, + ^{D^D^f'] 



At first order in ^ and second order in the B fields, this can be expressed as 



W 



iCAghr 



f^-'B^^U-^B^^ + 2{U-^)Bl^{U-^)Bl^{U-^)d^d, 



-AB^^AB;, - 2ABl,AB:^Ad^d^ 



where as usual □ = d^d^. 

We can write the renormalized expression of the first term of (5.12) as 

[A) = ^ J d^xd^y B;^{x)B-^{y)A^\n 

Whereas the second one is of the following form 

(B) = -2 j d^xd\d^u{dldlA^^)Bl^{x)Bl^{y)A,yAy 

= -2 1 d'xd'y Bl^{x)Bl^{y)A,ydldl J d\ A,„A,„ 

In order to evaluate this expression we must apply CDR in momentum space, arriving to 

1 



^yu 



Q, U L^xu^uy 



2™2 



4(47r2) 
-A{x-y) 



\n{x — y) m 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



(B) = 2 d^xd^yBl,{x)Bl^{y){Axydxd^A) 



(5.16) 
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Remembering the CDR renormalization of Adxd^A listed in (2.8) we can obtain the final renormalized 
expression for (B) as 



(B) = -^Jd'xd^yB%{x)B^^^{y)A'\n-j^Jd^xd^yB;,ix)Bt,{y)d,d,A 



(5.17) 



Adding up the two results 

{A) + {B) = -^^J d'xd'y Bl^Bl^d.d^A (5.18) 

This can be written at explicit second order in B fields as 

J B-^,{x)B'l^d^d^A = Jid^B-, - d,Bl){d,Bl - dxBDd^d^A + 0{B^) 

= - j B';Xx)B:{y){d^d.- 5 ^,U)U A + 0{B^) 

(5.19) 

Finally, remembering T = ~W + j Jcj) we can obtain the linear dependence in the gauge parameter ^ 
of the contribution at one loop in a generic gauge to the background two point function as 

Tc = -^^)Jd'xd% Bl{x)B'i{y){d,d. ~ S,,D){DA) (5.20) 

5.2 Two loop diagrams 

Now we follow with the two loop contribution to the background field selfenergy. Diagrams are those 
of figures 4 and 5 ((a) to (k)). First as an example, we will renormalize diagrams (a) and (k), as these 
two diagrams are an example of the two different types of divergences that we can find: Diagram (a) has 
nested divergences, whereas diagram (k) has overlapping divergences. The explicit evaluation of the rest 
of the diagrams is presented in appendix C. 

After discussing in detail those two diagrams, we will list the renormalized expressions of all of the 
two loop contributions that add up into the two loop renormalized background selfenergy. 

5.2.1 Diagram (a) 

This diagram has the following form 

< B^^{x)Bt{y) >, = _2g^faecfbcdfgdffgfe J ^4^^4^ ^^^(^ -d;){dy~ 5)A,„ X 

x(9^A™)A™(a^A,„) 

(5.21) 

We can rearrange this expression in terms of the integral I^ previously defined as 

< B;^{x)Bt{y) >a = -/Cl5"''[4a,(Aa^/i)-a^5,(A/i)-4Aa,A/'] (5.22) 
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v ftQi - ■ i.pQ .. 






Figure 4: Two loop YM diagrams a-f 



Now, in order to renormalize, we have to substitute these expressions with their renormahzed values 
(obtained in section 3.1), arriving to 



<Bl{x)Bl{Q)>aB, 



-Un^xHP 



32(47r2)3 

\\n^xHP + ^\nx^AP 

+(5^. □ □ ^ T^ 

x^ 



■Inx'^A'P 



(5.23) 



5.2.2 Diagram (k) 



In order to obtain all the contributions that form this diagram, the Mathematica package 'FeynCalc' 
was used, so that all the index contractions were performed by the computer, being its output the final 
relevant expressions that need to be renormahzed. The contributions shown here are those that have a 
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sAa; 







(j) 




Figure 5: Two loop YM diagrams g-k 



divergent part, omitting those terms that are finite. 



<i?;:(x)B^(y)>fc=icl<5"^ 



-16Sf.,dld:H[l,dxd, 



1,11 



20d:dlH[l,d^.dx; 1,1] 



~12Ad^dlH[l,dx; l,d^]+72dlH[l,d^dx; 1,9,.] 
+565^,dfd^H[l,dx ; l,d,] - 72df,H[l,d^D ; 1, 1] 
+20dp^H[l, □ ; 1, 1] - 144H[1, d,,D ; 1, d,] + 
+72d'^^H[l,D ; 1,5,] - 72d^H[l,d^dx ; 1,9a] + 
+Mdp-,H[l,dx ; l,dx] - 72H[l,d^dx ; 1,9,9^] 4 
+40ni7[l,9^9, ; l,l]+32niJ[l,9^; 1,9,] + 
+I6d^^,ni/[i, D ; 1,1]- l6Sf,,nH[l, dx ; l, 9a] ] 



(5.24) 



We now have to proceed as in the QED and SQED examples. Expressions that have a d'Alembertian 
can be put in terms of /^, and their renormalization is straightforward. The rest of the integrals can be 
found in the list of section 3.2 (or can be easily expressed in terms of integrals of that list). Finally , we 
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arrive to 



<Bl{x)Bl{Q)>uR = 






32(47r2)3 



f In'x^M^ + f Inx^M^ 
+<5^. □ □ ^ ^ 



(5.25) 



5.2.3 Two loop final results 

Proceeding in the same way as for diagrams (a) and (k) (see appendix C), we can obtain the renormal- 
ization of the rest of the diagrams, that we list here. 



(a) 
{b) 
(c) 
id) 
(e) 
(/) 
(5) 
(h) 

(i) 

(j) 
(fc) 



1 l„2^2»,f2 



2 A,f2 



1 1^2 ^2 A,r2 I 11 1„ ^2 A,r2 



v°°- 



■ In^ x^Af 



Inx^M^ 



- V°° 



f In^x^M^ + ninx^M^ 



a^a, □ ^ ^ 

-f In^x^AT^-flnx^M^ 
a^9.n^ ^ 

ilnx^A/^ 

(J^^nn^ + ... 

-Ih-ixHP 

'^M^°°^^2 + ■•■ 

i^^,^, - (5p.n)n 2 + ■ ■ ■ 

i In^ x^M^ - i Inx^M^ -i In^ x^AP - ^ hix^AP 
d^d^a^ ^ + S^^nn^ ^ 



a.9.n- 



^In'x^Af^ + flnx^A/^ 



5^,nn^2 



^In^x^A/^-^lnx^A/^ 
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I In^ x^M^ + f lnx2A/2 -| In^ x^A/^ + if \nxHP 
d^d.a^ -^ + S^^DD^ -^ + 

X'' x^ 

-^In^ xHP - i| Inx^M^ i In^ xHP + || Inx^M^ 



-^fi'^iy'- 



'--'fiiy '- 



1 1^2 ^2 ,^2 ^ 1 ,„ ^2 ,,,2 _ 1 lj^2 ^2^^2 _ 1 ln^2^2 

+ i5^.nn^ ^ + 



„■ In^ x^M^ + i lna;2Af2 -i In^ x'^AP - ^ Inx^M 
d^d^a^ ^ + S^^nn^ ^ 



x^ 

27 i„2 9 ,^2 45 l„ ^2 ,,r2 



a^^^.n " 



^In^x^AT^- ^Inx^Af 



-(5„^nn 



^Wx^M^ + f Inx^Af^ 



+ . 



(5.26) 



In this hst all the expressions have a common factor of 32(4^^)^ .9 "^^l'^"''- 

Finally, adding up all the terms we obtain the divergent part of the two loop renormalized expression 
for the background field selfenergy as 



4*^2 sab 



9'CIS' 



<B;{x)Bt{0)>n - -^(^;^v-M 



^" {d^d,~d^,D)D 



Inx^AP 



(5.27) 



5.3 Renormalization group equation 

With the previously obtained expressions for the one and two loop corrections of the background field 
propagator we can easily obtain the first two coefficients of the expansion of the beta function. 
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5.3.1 RG equation for quantum fields 



We have to consider the rcnormalization group equation for the quantum field propagator, as with this 
we can find the value of 7^, which is the function that takes care of the running of the gauge parameter. 
The only relevant coefficient of the expansion of 7^ in our two loop background calculation will be the 
first one. so that we only need the one loop correction of the quantum gauge field propagator. 
First, we write the effective action as 



d^xd^y Al{x)Al{y)T(^}-\x - y) + 0{A^) 



(5.28) 



If we consider the part of the Yang-Mills lagrangian which depends only in the quantum gauge fields, 
we can see that in a generic gauge is of the following form 



1 



1 



With this, the free action plus the gauge fixing term becomes 



Sq + Sg,f, 



1 



d^xd^y A^ {^5^,U5{x) - ^d^dj(x)) A^ + 0(^3) 



(5.29) 



(5.30) 



(2)ab 



Using the one loop calculation (5.9), r/tJ can be written as 



p(2)a6 



-5^,,u5{x) ~ £,d^,du5{x) + S'^^d^d, - (5^.n) 



5g^CA _\nx^M'^ 



-D- 



487r2(47r2) a;2 



-Six) 



727r2(47r2^ 
Inserting this into the RG equation 



0(/) 



.^^^+^(^)|+^4~'"". 



T^'^'"'k=o=0 



(5.31) 



(5.32) 



it is easy to obtain the desired result (along with the one loop coefficient for the anomalous dimension 
of gauge fields 7^) 



lA 



5Ca 2 
'24^^ 
5Ca 2 



487r' 



:9 



(5.33) 



5.3.2 RG equation for background fields 

If we define 

r(,2j-''(x) = (9^a.-5^,n)5-''r(2)(x) 

the equation we need to consider is 



^^alj+^(^)|+^4~'"". 



r(2'k=o = o 



(5.34) 



(5.35) 



If the background field is redefined as B' — gB this implies 75 = (due to the fact that in the 



background field approach the charge and background field renormalizations are related: Zg = Zg 



1/2 
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[15]). Then, with the one loop contribution (5.5), the gauge fixing renormahzation (5.20) and the two 
loop contribution (5.27) the effective action for the background fields is 

^ ' g2 "(■'=) + 48^2(4^2) ^ 3,2 + 72^2 ''^^^ + 3^2 "^^^ + 

+ -7 4^1=1 ^ h... 



2(47r2)a 



(5.36) 



Using (5.36) and the previously obtained expansion for 7^ (5.33) into the RG equation (5.35), we can 
straightforwardly evaluate the first two coefficients of the expansion of the beta function to be 

ncA 



/?: 



487r 



2 



17(^2 
These results agree with those previously obtained in the literature [15]. 

6 Conclusions 

In this paper we have shown how the results of CDR, found for the one loop case, can be very useful in two 
loop calculations. In those cases, the use of one loop CDR allows us to fix all the ambiguities that appear 
at one loop level, which implies that the coefficients of the divergent parts of the total expression are fixed 
a priori. No Ward identity is needed to be used to relate different scales. The procedure distinguishes 
between the cases of diagrams with nested divergences (where we renormalize from inside to outside the 
diagram, applying first CDR and after that normal DR) and diagrams with overlapping divergences. In 
order to deal with the latter, we have found a list of renormalized integrals that can be applied to various 
two loop two point function calculations. 

As an example of all this, we have found applying our procedure the two loop beta function coefficient 
of QED, SuperQED and Yang-Mills theory. 

Acknowledgments 

I would like to thank J. Mas for useful discussions and critical reading of the manuscript. Also I'm very 
grateful to M. Perez- Victoria and J.I. Latorrc for sharing their insight with me on CDR and the two loop 
differential renormahzation of QED, respectively. Finally, I would like to thank M.Gomez-Rcino for her 
critical reading of some parts of the manuscript. 

A Integrals with overlapping divergences 

Here we present the explicit calculation of the list of integrals of section 3.2. As there, we define the final 
results in terms of a variable z = x — y and denote 9^ = df^. Also, in each final expression . . . stands for 
the finite contribution that we are not taking into account. 
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• H[l,l; 1,1] 

i^[a^,i; 1,1] 
H[l,dx; 1,9a] 

These integrals are obtained by means of Gegenbauer Polynomials [1, 5]. 

• dlH[l,d^; 1,9a] 

Contracting (D.l) with S,yx wc find 

d^Hll, d^ ; 1, 5a] = \dlH[l, 1 ; □, 1] - i/[l, d^ ; 1, D] - ini7[l, 9^ ; 1, 1] 
R 1 ilnz^Af^ 

a^n^ . + ... (A.l) 



32(47r2)3 ^ 

H[d^dx,dx- 1,1] 

In this case, no integral relation is used 



\dldlH[l^dx- l,l]-i 

ia,(A/i)-^(a,<5) 

1 -iln^z2A/2_lnz2jV/2 



a;^ajif[l,aA; l,l]--a^i7[l,n; 1,1] 



d^u^ + ... (A.2) 



32(47r2)3"'^ 

dlH[l,dx; d^d,,l] 

H[l,dx] 5a 5^,1] 

First of all, the third integral (3.10) will be evaluated with relation (D.l) 

H[l, dx ; dxd,, 1] - ^d^xHih dx ; 1, d,] + ^dlH[l, 1 ; 1, ^a^,] + ^d^Hil, dx ; 1, ^a] + 

+ ^d^dlH[l,dx; 1,1] (A.3) 



Using the previous results 



1 iln2z2M2-|lnz2^/2 



H^[l,dx;dxd.,l] = ^^(^^9,n^ ,2' +■■• (A.4) 

But this integral along with the other two, can be obtained with other method. We can apply the 
CDR decomposition into trace part, traceless part and local term (2.7) to the divergent subdiagram 
{d^d^ Ayu) Ayv Auv Let us consider this in the general integral 



/ 



d-'ud^v A„,(9^A,„)(5|'a^A^„)A^,A„„ = 

- '-S,^iAd,I^),-^d,Al-^I,,.n (A.5) 
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where Ip^a stands for the traceless part. The one loop ambiguity fixed by CDR is reflected in the 
second term of (A. 5) that at two loops has becamed a logarithm of the scale. In the renormalization 
of the traceless part normal differential renormalization will be used, leaving ambiguities (local 
terms) not fixed. 

The expression for Ipi^cr is 



Wn = B^-^-\a"^M + ^A-\b) 






(A.6) 



or in the form of the integrals being discussed 

3 



'AA/j R 



{A7r'ri3A-B)d,Al 









(A.7) 
(A.8) 

(A.9) 



The value of (3^ — B) is easily obtained using (A.7), because this corresponds to integral (3.10) 
that was obtained previously. 

I.e. 



d'^ud'^v A„,{dlA,,){dldyAy^)Ay,A^ 



32(47r2) 
1 



1 iln2z2M2-Ilnz2Af2 

— a„n^ -^ 



3"M" 



32(47r2) 



:d,U^ 



Hn'zHd^ + ^lnz^AP 
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(47r2)£ 



ipea R, 



(A.IO) 



which implies that 



i^A-B) ^ ^ 



(A.ll) 



With this result the evaluation of (3.12) and (3.14) are straightforward 

1 



9fH^[l,9,, ; d,dx,l] 



^dAAd.n,- 256.2 



dpduAfi — /, ^^^ dxl^xu . 



(47r 



27^ClX-lf_iXu R 



1 



32(47r- 



2^3 



dpd,a» -^ + A nn^ 



(A.12) 



dfR^'ll^dx; 9^9., 1] 



-]s,^dx{Adxl')R - 4^0Al - j^dxh,. R 



2567r2 



1 



32(47r2)3 



-ilnz2Af2 

dud^n^: — + Sn^aa 



(47r2)5 

1 In' zHP 



llnzHP 



(A.13) 
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. dlH[l,dx; hd^d,] 



Using integral relation D.l 

dlH[l, dx ; 1, d^d,] = dfH[l, dx ; d^d,, 1] - dld^^H[l, dx ; 1, d,] - d^,dlH[l, d^ ; 1, d^] 
- d,aH[l,d,r, 1,1] 



32(47r2)3 



dud.a^ 



^Inz^AP 



1 1„2 ,2A,f2 I 3 l„ ,2»/r2 



,5^,nnA 



In-'z^M^ + llnz^M^ 



(A.14) 



Considering (3.12) and applying (D.l) 



d-xH[l,d^; d,dx,l] 



= ini/[i, d^ ; 1, aj + ld-^diH[i, 1 ■ 1, d^dx] 

+ Id^dlHll, d^ ; 1, dx] + Id^OHil, 8^ ; 1, 1] 



(A.15) 



Remembering previous results 



1 



1 



dlH^[l, d^ ■ d.dx, 1] = ^aH^[l,d^ ; 1, a.] ^ 33(4^2)3 -M^-- 



ra^a^D- 



iln2z2A/2 + ilnz2M2 



So that 



nij«[l,9^; 1,9,] 



32(47r2) 



1 , -^Inz^AP 
5^,nn^ 



(A.16) 



(A.17) 



Using (3.14), (3.15) and the identity 

DH[1, 1 ; d^d,, 1] = a^ff [1, dx ; d^d,, 1] + d^xH[dx, 1 ; 9^9,, 1] 



(A.18) 



we can easily arrive to 



„, , 1 iln2z2A/2 + |lnz2M2 



32(47r2) 



3 ^'^ 



(A.19) 



• 5^JI[1,1; dxd,.,d^] 



Using (D.3) and (3.16) 



dfH[l,l; dxd,,d,,] = ^□i/[l,l; a^9.,l]-ia;i9^a^i7[l,l; dxA]-\dyaH[l,l; d^,l] 



1 1, 2 ,2»,r2 I 3 1 ,2 7i,r2 



/J 1 iln z^A'/^ + llnz^M^ 



32(47r2)3 '"' 



(A.20) 
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. ajF[l,l; d^d,,dx] 



With (D.3) and (3.9) 



dlH[l, 1 ; d^d^.dx] = -diUH[l, 1 ; 9^, 1] + -^f a^a^i/[l, 1 ; 9a, 1] + T;dldlH[l, 1 ; dxd,, 1] 



32(47r2) 



1 \\n^z^AP + l\nz^AP 
dnd,U& ^ 



2^3^^^!^ 



(A.21) 



In this case the CDR decomposition into trace part, traceless part and fixed local term (2.7) will 
be used again, as in (3.12) and (3.14) 



H[l,d^,d\ ; dyd\,l] 



-^{^^^^,I^)R - \[^{d^d, - iv°)^% 



1 



r(Aa,,9,A)j 



1287r2 ^""''^ 

64 

J^^l^,\u\ R 



1287r' 



^[A(9^9,-iv°)^]i? 



(A.22) 



where /^ai/A stands for the integral with the traceless parts. This was calculated in [2], and the 
result was found there to be 



64 



(47r2)5 

Adding up all the terms, it is easy to arrive to 

1 



^/pA^A B. = 96(4^2) ^A'^'^^fl + 7^'^M^^fl 



13 

48' 



ii«[i,a^aA; d^dxA] 



32(47r- 



2^3 



z^ 



-^In^zHP-^lnzHP 

V°°^^ r^ 



In this case, applying (D.l), (3.9), (3.10), (3.12) and (3.19) we get 



(A.23) 



(A.24) 



H[l,df,dx ; 1,i9a9,. 



H[l, d^dx ; d,dx, 1] - dlH[l, d^dx ; 1, d,] - dlH[l, d^dx ; 1, 9a] 

d:diH[i,d^dx-, 1,1] 



1 



32(47r2)3 



1 1„2 ^2 7i/r2 _L 49 l„ ^2 7, ,r2 



a^a.n^ 



In^z2^j2_^4yij^^2^,j 
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-+ 



5,,,un^± 



1 ln'^2^2_ llij^^2^j2 



(A.25) 
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B Expansion of n^Ti^ 

The complete expanded expression of the two loop diagram b of QED is 



n(2&) 



-86^,aH[l, dx ; dx, 1] + l6d^^H[l, d, ; □, 1] - 86^,dlH[l, dx ;□,!]- 
~16d^^H[l, dx ; dxd,., 1] + 16dfH[l, dx ; d,,d,, 1] - mdfH[l, d^ ; Oxd,, 1] 
~16df,H[l, D ; a„ 1] + 85^,dlH[l, n ; dx, 1] + 16dlH[l, dxd^ ; 9,, 1] - 
-l&dlH[l, d^d, ; 5a, 1] + 16a^iJ[l, a^a^ ; dx, l] - 16i7[l, □ ; d^d,, 1] + 
+85^,i7[l, □;□,!]+ 32i/[l, a^^A ; a.^A, 1] - 16ff [1, d^d, ; □, 1] ] 



(B.l) 



C Two loop diagrams of Yang-Mills theory 

We remind that in these diagrams I?^ denotes a derivative acting over one of the external fields, and 
again . . . stand for the two loop finite contributions that we are not taking into account. 



C.l Diagram (b) 

C.l.l Definitions 

When obtaining this diagram some quantities are defined in order to simplify the expressions (as we 
defined I-^). They are 



^2 

UV 



/O(x-y) = j d^ud% A,uAy^{A^^d;Auv) 



(C.l) 
(C.2) 

We will detail the renormalization of the terms formed with these integrals that appear in diagram 
(b): AD/o, Adf,dJ° and A/o^. 

• /" has UV and IR divergences. This expression was studied in [4] and evaluated as 



7^° 



^^^^ [In' xHlfj, + 2lnxHlfj,{l - InxHl^) + bji,] 



It is clear also that D/*^ = — I^, and it is not difficult to prove that 

1 



(^^'^^^^°)« = 32(4^ 



dnd^ 5 + Sn^n* -^ 



(C.4) 



(C.5) 



The renormalization of /" is straightforward, once we recall that CDR imposes /° ^ ~ h'^nl] 



2^tJ.^R 



With /|J^, applying CDR to the subdivergency we find 

II. = ^ I d^ud\, A^^Ay^{d,d^ ^ y,^D){Al,) 



2887r 



2 , d'^ud'^v AxuAyy{d^di, - Sf^^O)S{u ~~ v) 
1 



= ia^a./o - i<5„„n/"„ 



12^'^"^ ^^72(4712^"'^"'' 



d„d,. I d'^u A^u^yu + -,2(7 2] 
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or, which is the same 



(AV)fl 32(4^2)3 



C.1.2 Evaluation of diagram (b) 

This diagram is of the form 



d^d,^ 



Inx^AP 



V°- 



■lna;2M2 



x^ 



(C.6) 



< B;{x)Btiy) >, 



2 race rbed I i4 t4 a 

g J J I a ua V /\^. 



^2S,,D;+S,,{df,~d-J + 26,pDl 



xA'ii{u^v)A,y 



-25,„Dl + S^xid^- dy) + 2S,,xDy 



A, 



(C.7) 



Where A^i^ is the one loop correction of the quantum gauge field propagator, which we have found 
to have the following non renormalized expression 



A;t = g^cls""" [a^a,A2 - j^^da^ + 4a^(Aa,A) - 2<5^,a^(A5AA)- 



(C.8) 



It has to be noted that, in contrast with dimensional regularization, the renormalized one loop expres- 
sion can not be used here. The reason is that the indices of the one loop insertion the will be contracted 
in a future step, and one of the rules of CDR is to make first all the index contractions before performing 
the renormalization. So that, we are only allowed to insert the bare one loop propagator. Then, the 
diagram can be put as 



< Bl{x)Bl{y) >t= -g^CAS""'! d^udS 



-49;^ a^ (A^„Ap^(ii - v)A^yA^y) 

+A6f,^d^d^ {A^uApa{u - v)AyyA^y) 

+AUd; ~dl)Ap,{u - v)A,y{dl~ 5)A., 
+4n (A^„A^^(m - w)A,„j,A^j,) 
-4.dldl {A.^uAf,„{u - v)A^yA^y) . 

Straightforwardly we write this expression in terms of the previously defined integrals as 

< Bl{x)Bl{Q) >bR = g^ClS"- [^2Ad,,d,{Ad,d,l") + Udpd,{Aai") + 32a^a,(A4"J + 
+ 125^,d^dp{Adpd^l'') - 16<5^,n(An/") - 165^,9^9, (A/;,) + 
+ 205^(A9„n/") - 20Adpd,Dl" + l2n(Aa^a,/") - l6n(A/;;j]^ 

And with the renormalized forms of /^ and /" is easy to arrive to the final renormalized expression 



<B;ix)BtiO)>bR 



g'^ClS"'' 
32(47r2)3 



d.d^a^ -^ 



A nn-2 2 



(C.9) 
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C.2 Diagram (c) 

The expression for this diagram is 

<Bl{x)Bl{Q)>,R 



4 ^acx ^xed rbdy syecx A 3 



„4^^2 rab 



g-'C^AS^'S^^A' 



g'cls-" 



32(47r2)3 f"" 

C.3 Diagram (d) 

This diagram is similar to the previous one, and its result is 

9 



<5„„nn^ 



(CIO) 



<Bf,ix)Bt{0)>dB 



C.4 Diagram (e) 

The bare expression for this diagram is 



g^C'lS'^'S^^A^ 



2 A,f2 



32(4712)3'^'"'^^ 



(C.ll) 



< B^,(x)Bt{y) >e 






+ 25,xDy] Ayu 

= ^Qg^ClS'^^dldl - 5,.U) j d*u Al^Al 

The renormalized expression of the integral is easily obtained as 

1 1 



26,,Di + s,^{dy-dy)+ 



I 



d\ 



{x — u)"* u 



4„4 



7r2 In^ a;2M2 

— □ ^ 



So that 



■i 1 2 2jir2 

<Bl{x)Bt{Q)>eR = ^7r^T^g'cl6''''{d^d.-S^,a)D^^ — 



8(4^2) 



(C.12) 



C.5 Diagram (f) 

This diagram is 

< B",{x)Bl{y) >f 



^gA^acj.dfjdcejbeJ ^ ^4^ A^ (9;: A„,) (9^- d^A^y + 
^g4fafcjecdfbf.f.ed /'rf4^^^^(^ -9^ A,, (9," A,„) A„, 
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Operating 



<Bl{x)Bl{Q)>jR 






32(47r2)3 






Sa.^n^ -^ 



+ ... 



(C.13) 



C.6 Diagram (g) 



< Bl{x)Bl{y) >, 



It is easy to arrive to 



_2gifac.f.fdfecdjbfe l^i^ 5,,, A„,(9^ A,„) A„^ [-25,Ai?^ 



6x4dy-dy) + 2S.,Di 



Axu 



<B-{x)Bt{0)>,R = -g^ClS'^' 

32(47r2)3 



-9^(Aa,/i) - A9^9,/i - 5^,0,(AaA/i) 



5^9.n 



^ln'x2A/2 + i|lnx2A^2 



+ 



-^ln2a;2M2-Xina;2Af2 
(5,,, D D -^ -J 



(C.14) 



C.7 Diagram (h) 



< S^(x)i?^(y) >,.= 






X AW/^'^ 
X Ai;)/^/^ 



<5a.(9," - a;') + <5Ap(a;' - 9^) + s.^idi di) 



-2S,xDy + 5,x{dl- dl) + 2S,,Dl 



In this expression we use the convention of A^*^ 9^ A^-') = (9^A'^'))A(-'^. Evaluating aU the index 
contractions 



<B^^ix)BtiO)>,,n = -g'C^S^' 



g^C^S"'' 
32(47r2)3 



4^ 

-d^iAd^l') ~ 27A9^9,/i - 95^,dxiAdxl') 

Un'x^N'P+'^lnxHP 
d^d. □ ^ 2-^ + 



J^^DD 



-|ln2a;2Af2+ illna;2A/2 

4 4 



(C.15) 
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C.8 Diagram (i) 



This diagram and the next one have overlapping divergences as discussed with diagram (k). We wiU 
proceed as in that diagram, making use of the overlapping integrals listed in section 3.2. 
The bare expression for this diagram is 



and operating this becomes 

< B;{x)Bt{y) >,= -^g'ClS^' [ d;dyH[l,dx ; 9a, 1] - 2d-^H[l,dx ; dxd,, 1] - 

-2dyH[l, dxd^ ; d^, 1] + 4H[1, 0^0^ ; d^d^, 1] 

Finally, with (3.7), (3.10), and (3.20) the renormalized expression is 



<B"^ix)BliO)>,n 



32(47r2)3 



dndu □ -^ TT^ 



^WxHP + ^lnxHP 
5^.UU^ -^ 



C.9 Diagram (j) 

This diagram is of the following form 

< B';{x)Bl{y) >, = 2g^rfr''f''^jf^y j d\id^v A 

+ 25,pDl] {d;A^y){dy- 9"^)A,„(a:A„„)A 
Evaluating the index contractions this becomes 
< B:{x)Bt{y) >,= -g^C^S'^" [ 



-2d,.D; + Sp.{df,~d;)+ 



~Ad:dlH[l, d^dx ; 1, 1] - ^dldlH[l, dx ; 1, 9^] 
-4a^i7[l, d^dx ; 1, d,] - Ad^xH[l, dx ; 1, d^d,] ~ 
-Ad^,H[l,d^dx ; l,dx] + dp^H[l,dx ; 1,9a] - 
-AH[l,d^dx; l,9,9A] + 4niJ[l,a^; 1,9,] + 
+4nif [1, 1 ; 1, d^d,] - 2d^H[l, d^ ; 1, D] + 
+d;d':H[i, 1 ; 1, n] - 4H[1, d,a ; 1, d^] - 
-2d;H[l,d,a; 1,1]] 



Proceeding in the same way as in diagram (k) we find the renormalized expression to be 
< B-{x)Bt{y) >, 



gHJl6^ 
32(47r2)3 



■ ^ iln^x^M^ + ilnx^M^ 



-iln^x^M^-ilnx^M^ 
A nn ± ± 



+ 



(C.16) 



(C.17) 



(C.18) 



(C.19) 



(C.20) 



(C.21) 
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D Integral relations 

To deal with overlapping divergences, the following exact integral relations are very useful 

dlj d\id^v A,„(9^A,^„)A,„(a^A,„)A„„ + 
+ Jd^ud\, A^u{dlA,,)Ayu{dldlAy,)Au^, + 

+di I d\d\ A^„(a;^A„,)A.y„(a^Aj,„)A™ + 

+di;dl / d'^ud'^v A^u{dlAxv)AyuAyyAuv 

dl j d'^ud'^v A^uAxv{dldlAyy)Ay^Auv 

= 2^1 j d^ud^v A„Axv{OAyu)AyyAuv - 
- / d^ud^v A:^uAxv{^Ayu)AyuAuv + 
+\dldl j d^ud\ A^,uA,^{dlAyu)Ay,Au. 



(D.l) 



(D.2) 



dl / d^ud^v A„A^,{dldyAuy){dlAy,)Au, 

\dldldl j d\d\ A,uA,^{dlAuy)Ay,Au, 
+ldlu j d^ud\ A,^A,,{dlAuy)Ay,Au. 
^in j d^ud\ A,uA,^{dydyAyu)Ay,Au, 



(D.3) 
To prove these relations, wc have only to perform the derivatives and expand the different terms. 
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